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SUMMARY 



A series of calculations was made to determine the probable 
dynamic wing street of two large airplanes in atmospheric gusts. 
These calculations were undertaken oonourrently with a more general 
investigation* still incomplete, from which it appears at the 
present time that the calculated stress agrees well with measurements 
on flexible wing models in the gust tunnel. 

The results of these special calculations indicate that in 
both isolated and repeated gusts of probable occurrence the dynamic 
over stress is about 10 percent when referred to the present static 
design standard. 



INTRODUCTION 



The possibility of dynamic over stress in airplane wings upon 
encountering atmospheric gusts has been the subject of a number of 
investigations. The formulas and methods resulting Aram these 
investigations when applied to specif io cases as a rule showed 
overstressj whereas preliminary tests of a flexible-wing model in 
the gust tunnel consistently showed understress. In view of this 
oonf Hot it was felt that a further development of the theory and 
more comprehensive tests should be undertaken. - 

This work has been carried out to an advanced stage , and 
reasonably good agreement has been found between the theoretical 
and experimental results so far as this part of the work has been 
carried. The analysis of actual oases has, however , been largely 
confined to obsolescent designs in single gusts, in which oases 
understress rather than overstress was the outstanding result. In 
view of the rapidly changing trend in design; it was felt desirable 
to apply the results of the investigation to two airplanes currently 
in the design stage of development, and to extend the analysis to 
the effect of repeated gusts* This special analysis disclosed 
possibility of serious overetresei accordingly, the results are 
presented herewith for the information of those concerned* A more 
comprehensive report on the complete investigation is intended 
to follow* 



2 



For the purpose of this report, the airplane designs in question 
will be designated model A and model B. 



METHOD 



Inasmuch as a more detailed di-eeussion of the theory 1 will 
be presented in a later report, only a brief outline will be given 
here* 

An airplane flying through the air may be oonsidered in tha 
spanwise direction as a beam of nonuniform oross section supported 
on a yielding foundation. Slnoe a rigid solution of the problem 
of dynamic stress for suoh a beam is impractical, as pointed out 
by El saner in reference 1, an equivalent wing and spring system for 
an airplane was assumed as indicated in figure 1. The motion of the 
upper wing is adjusted to be that of the wing tip motion of the 
original airplane and the motion of the lower wing and fuselage is 
the motion of the fuselage of the original airplane under shear 
load fir am its wing. The equations of motion for this system under 
the influence of a single gust of the type shown In figure 2(a) 
are as follows* 



M^D 3 ^ + j^B^ ♦ KdW - 8 f ) = (D 



Itf D 3 ^ + \ £ Dfi f - KCV - 6 f ) = te" bt (2) 
e e o 



where the subscript w refers to wing quantities and the subscript f 
refers to fuselage quantities, and 



t 


time 




equivalent mass of wing 


"f. 


equivalent mass of fuselage 


D 


differential operator 


«ir 


apace position of wing (fig. 1) 


8 f 


space position of fuselage (fig. 1) 
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X 



damping coefficient = + Xf = (effective damping 
factor) mf-SV e 8 



m 



slope of lift curve, per radian 



P 



air density 
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wing area 



V 



airspeed 



The equivalent mass of the wing, U_ , is the mass which, if 

a 

placed at the wing tip of a weightless beam, would give the same 
deflection under unit acceleration as the distributed mass of the 



In the absence of knowledge of the actual deflection curve of 

a wing, an assumption must be made for the purpose of computing the 

damping coefficient* If, as will be assumed here, the deflection 

at any point is proportional to the square of the distance frcm the 

root, the damping of the wing tip motion is equal to the damping of 

the whole wing times the ratio of the deflection of the mean vertical 

velocity position to the tip deflection, or \ w = Th* effective 

e 

damping factor for the range of wing frequencies of interest is 
taken as 0*73 frcm the results of a preliminary analytical and 
experimental investigation. 

The total air load on the wing may be considered as consisting 
of two components, one of whioh results in bending deflection at 
the tip and the other of whioh results in a shear force at the 
wing root* The division of the air load into these components is 
accomplished in much the same manner as the division of mass of the 
wing* The shape of the air load on the wing was assumed in the 
oases of both airplanes investigated to be similar to that on an 



wings In this oase it is 




as shewn in figure 1. 



The equivalent mass of the fuselage, M f , is equal to the mass 
of the fuselage plus the actual mass of the wing minus • 



averags tapered wing end ijb waa divided so that A^te applied 
at the wing tip gave the same deflection aa the total distributed 
air load* For the oases under Investigation It oan be shown 
that 4_ = 0#25A where 

6 

A s at bt = 1 (3) 

te"" 

W total weight of the airplane 

1 -bt 
b = f at bt s 1 or the maximum value of the function te 

Since A^ and A^ appear as multiplying faotors for the 

equations in their solved form, An may be any arbitrary load 
factor increment i An is the load factor increment that the 
airplane would experience if it had no vertical motion as it traveled 
through the gust. 

The function te was taken as most closely representing 
the shape of the time history of air load on the wings in a gust 
as Indicated by tests in the gust tunnel. 

Equations (1) and (2) were put in operational form, combined, 
and solved to give the following results t 



6 f s e"*** (ox oos R 3 t + o a sin Bat) + ^e^* 3 * + o 4 

♦ Kite'* 5 + fce'** (U) 



*w = ~-D 3 6f * ~-Wf * «f - ~- **~ bt (5) 



An expression for the normal acceleration of the airplane when 
the wings are held rigid may be derived \yy letting (V • 8 f) 
equal eero, so that D 3 6f * D?^ and DS f s DC^ in equations (1) 
ax& (2). With this restriction equations (1) and (2) may be oombined 
to becomes 



5 



where D*^ s in,,, the aooeleratlon of the rigid airplane. 

The gust gradient distance Hi ie determined by the time 
interval from the start to the peak of the br^ curve. Oust- 
tunnel teste of rigid wing models substantiate this method since 
the normal aooeleratlon on the model reaches its maximum value at 
the same time the gust velocity of a linear gust gradient reaches 
its maximum value. 

For the purpose of these calculations the wing deflection is 
assumed to be proportional to the wing stress* The dynamic stress 
may therefore be given in torms of the ratio of the dynamic wing 
tip deflection (&w - 6f)» which is designated 64, to the static 
deflection • 6 fl ^ under the same ^condition of load. In this oase fi at 
may be determined from the equation 

Any x 0.25(W tot ) An, 1L 
«st = , 2 T— 1 (?) 



The time history of reactions calculated for a single gust are 
basic curves to which the principle of superposition may be applied 
to determine the reactions for repeated gusts. This method is 
substantiated in reference 2 where an analogous problem in eleotrioity 
is presented. 

in alternative method of determining the reactions for a 
repeated gust would be to apply a suitable repeating force function 
to the equation. This method would be tedious, however, since each 
gust combination would have to be calculated separately) whereas, . 
by superposition, three basic sets of curves representing different 
probable gust gradient distances may be used to give an almost 
unlimited number of combinations of two gusts. 

There is no direct mathematical connection between the gust 
velocity distribution and the forcing function, since the forcing 
function was determined by inspection of gust-tunnel records. The 
distribution must, however, approximate that shown in figure 2(a). 
The gust velocity distribution for a repeated gust therefore may 
be obtained by adding together two gusts of this type as shown , in 
figure 2(b). The result of this addition may be of the type shown 
in figure 2(o) or 2(d). 



Although theeffebt of airplane stability is not taken into 
account directly, it is felt that the form of the forcing function 
la same degree takes into account this effect. However, it is 
believed advisable, because of the limitations introduced by lack 
of direct consideration of the pitching motion, to limit the number 
of superpositions to tiro gusts* 



CALCULATIONS AND RESULTS 



The conditions and basic constants used for the calculation 
of the reactions of model A and model B are given in table I. 

In the practical application of dynamic stress calculations, 
it is essential that the true nature of the force causing the 
stress be known* In the present instance, this is equivalent to 
saying that the gradient distance of the gust be known, since it is 
this distance that determines the nature of the forcing function* 
Fortunately, it can be stated that the most probable gradient 
distance associated with tho largest gust accelerations is about 
10 chord lengths and that the severity of the acceleration falls off 
rapidly with a departure of the gradient dlstanoe from this value. 
This fact has been well demonstrated by extensive flight tests on 
airplanes varying in size from about 1500 pounds to 65,000 pounds, 
and can be demonstrated by physical considerations* In order to 
obtain results of practical interest, therefore, the calculations 
for models A and B have been carried out for gusts having a gradient 
distance of 10 chord lengths* In addition to these calculations, 
however, further calculations have been made for gradient distanoes 
of about h chord lengths and of about 20 chord lengths in order to 
extend the range of the analysis. 

Time history curves of the reactions of models A and B to 
single gusts having gradient distanoes of about h* 10, and 20 chord 
lengths are given in figures 3 to 5* 

Sample time history curves for a repeated gust composed of 
two equal and opposite 10 chord- length gusts arranged relative to 
ono another to give maximum negative wing deflection are given in 
figures 6 through 8. 

Figures 9 through 11 show for single gusts the variation with 
gradient dlstanoe of the dynamic stress ratio or ratio of maximum 
dynamlo wing deflection to the static deflection. They also show 
the variation of the ratios of maximum wing-tip and fuselage 
accelerations to the maximum accelerations for the case of the' rigid 
wing. 
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- — Since tha interest is primarily in dyxaunio . stress. rather than 
in accelerations, the superposition of the calculated basic curves 
was done with a view to determining the maximum over stress from the 
combination of the reactions of two gusts. It was found that maximum 
values of dynamic stress ratio occurred when the repeat gust was 
of a negative sense in relation to the first gust and that the 
sequence period or distance H 3 (see fig. 2) had a pronounced 
effect on the result* The maximum value of. A^, occurring in the 

whole sequence was used to determine the static defleotion. Table II 
presents selected oases which were the most serious of a number of 
combinations examined. In this table the results refer to repeat 
gusts having the same values of An as the initial gusts. It should 
be noted that the true gust velocities of the Initial and repeat 
gusts for any condition are therefore the same, and independent of 
the value of E since, it will be remembered. An is the load 
factor increment that the airplane would experience if it had no 
vertical motion as it traveled through the gust. 



DISCUSSION OF RESULTS 



The results of the calculations for model A in a single gust 
(figs. 9 and 10) and for model B (fig. 11) show that the dynamic 
stress in all oases increases tram, an under stress for a gust of 
long gradient distance to moderate to high over stress for a short 
gradient distance. The general shape of the curves seems to indicate 
that even higher over stress would occur in sharper gusts than those 
examined* However, the lag in development of lift in sharper gusts 
would preclude suoh a result f since even for an infinitely sharp 
gust, the forcing function would have a character similar to that 
for h chord lengths. 

Differences in the values of dynamic stress between the two 
conditions of model A would indicate that the change in forward 
speed and weight condition of the airplane has a pronounced effect. 
The more general analysis previously. mentioned has shown that, in 
general, the dynamic stress increases a small amount with Increase 
in velocity but that the effect is small. The change In weight and 
weight distribution oroduoes a substantial ohange in dynamic stress 
and, In general, a change in these qualities of such nature as to 
reduce the wing frequency' will result in an Increase In the value of 
the dynamic stress. 

"When the basic curves in thifc paper cure superposed to obtain 
the maximum overstress from their combination, it is seen from the 
results presented in table II that there is no definite correlation 
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between the effeot of gradient distance of the first and seoctad gusts 
and the distance between them H3. This lack of correlation results 
from the influence of ■ certain other factors, such as the relation 
between the time to peak acceleration and the period of wing 
vibration, which complicate the problem when the reactions to one 
gust are superposed on those to another gust. 

Examination of the values in table II shows that substantial, 
overstress exists for all the combinations of gusts presented and 
that the addition of a short gradient gust produces the largest 
value* As indicated in the preceding disoussion, the value of 
overstress becomes large in all oases as the gradient distanoe of 
the gust is decreased. However, before an estimate can be made as 
to whether the overstress on the airplanes in question will be 
serious, it will be neoessary to consider the effeot of the intensity 
and size of gusts and their spatial distribution in the atmosphere. 



APPLICATION OF RESULTS 



It should be emphasised that any process chosen at this time 
to appraise the significance of the dynamic stress calculations 
with respeot to practical design questions must be viewed with 
considerable suspicion. The subject is entirely too broad and 
involved to be given adequate treatment here although, in order to 
preclude improper application of the results, a brief disoussion 
seems necessary at this point. 

The present design oritiericais for gust loads are based on 
thousands of hours of acceleration-air speed data obtained with 
Y-Gt recorders installed on transport airplanes. The analysis 
of these data has indicated that a reasonable value of the effective 
gust velocity to be used for design purposes is 30 feet per second. 
As previously mentioned, the most probable gradient distance 
associated with this gust velocity has been found from separate 
investigations, such as that reported in reference 3, to be about 
10 chord lengths. As shown further in these investigations, the 
gust velooity measured on a given airplane increases from a negligible 
value to a maximum as the gradient distance is increased from 0 
to 10 chord lengths, but with further increase in gradient distanoe 
the measured gust velooity tends to fall off. 

In view of these facts it is felt that the overstress indicated 
in the short single gusts may be disregarded for the reason that 
the total Btress will be less than the stress for the more important 
gusts even without overstress. 
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It. is- ftirther; feit, howevexv ±hafc -the: amount of over ■tread 
shown for the airplanes In question in single gusts when the 
gradient distanoe is 10 chord lengths (figs. 10. and 11) should 
he added to the design stress oaloulated on the usual basis. 

, So far as repeated. gusts are concerned, there are no data 
available concerning the size. Intensity, and probability of 
occurrence of such gusts, einoe no practical method of analysis 
of available reoords to obtain such information has been found to a 
date* It is necessary, therefore, to appraise the probable over- 
stress in repeated gusts on the basis of measurements of the . 
structure and the frequency of single gusts* For this purpose let 
it be assumed, for the sake of illustration, that any airplane will, 
during its life, encounter only one single gust of the else and 
Intensity corresponding to the present design effective gust* With 
this assumption, It follows that each gust of any probable combination 
of Repeated gusts must have less intensity than the single design 
gust* Such a restrictive assumption is not really necessary, as 
the same argument applies in a relative or qualitative sense If the 
design effective gust is encountered several times during the 
airplane life* 

In order to determine the reduced value of the repeat gusts, 
recourse Was had to gust measurements made on a number of airplanes. 
Reference 3 presents data of the type used for this purpose* These 
data were applied on the assumption that the maximum Intensity of 
the three gusts next In Intensity to the maximum gust appearing in 
the data would bear approximately the same ratio to the intensity 
of this maximum gust as the intensities of the probable repeated 
gusts would bear to the single design gust. With this assumption 
the ratios were found from the several sets of data to be 0.75, 
0*66, and 0*6l* Thus the intensity of the individual gusts in a 
repeat combination may, depending an the data used, be 



TJ rpt " 0.75 V et0 * 

Referring to figure 2, this means that the repeat gust combination 
would appear as In (o) with the value of equal to 0.75 

(eto*) times the value of IT, appearing In (a). 
■ ■ . *max. 

The ratio of the maximum stress experienced in the repeat 
gust sequence to the static stress usually oaloulated for the standard 
design gust may thus be obtained by multiplying the ratio of the 
m a xi m u m dynamic stress in the sequence to the static stress In the 
first phase of the repeat combination the ratio of U r pj. 
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to T^^. Mote that this process differs, fs-om that employed in 
obtaining the dynamic stress ratios of table II in which the 
naxlmum value of A zip in the whole sequence was used to determine 
the static deflection. 

Confining the remainder of the analysis to repeat gusts having 
values of H a and • E a of 10 ohord lengths, but utilising any 
value of H 3 . oases may be selected out of table II for further 
examination. The oases se looted are those marked with an asterisk 
in the first column of table It. The maximum dynamic stress in these 
oases, referred to the static stress in the first phase, are as 
follows I 

Model Condition 8 d/ 6 st 

A I . 1*58 

A II 1.6U 

B I 1.62' 



These values, multiplied by the repeat-gust Intensity ratios 0«75t 
0.66, and 0.61, previously explained, have been plotted in figure 12 
against the intensity ratio. On the whole, figure 12 indicates that 
8 erne dynamic overstress is to be expected in repeat gusts relative to 
the present design standard* Taking a mean value, the overstress 
is about 10 percent, or about the same amount found for the single 
gust. In view of these results, it is felt that the design strength 
for these airplanes in the gust condition should be increased 
about 10 percent. 



CONCUJSiaSS 



The analysis of the results of the calculations for models A 
and B indicate si 

1. The dynamio stress in the airplane wings from encountering 
gusts in the atmosphere increases as gradient dlstanoe decreases. 

2. For the established important gradient distance of 10 ohord 
lengths, the overstress in a single gust is about 10 percent. 
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- 3* Hhmty probable types of repeated. gutft. are, enoountered* the. 
airplane wings may be stressed about 0 to 20 percent beyond the 
stress allowed by the present design gust load factor. 



iA&gley Memorial Aeronautical laboratory, 

National Advisory Committee for Aeronautics, 
Iangley Field, 7a. 
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TABLE I.- CONDITIONS AND BASIC CONSTANTS OF MODEL A AND MODEL B 

Model A 





Condition I 


Condition II 


Weight, lb 


6e,soo 


102,000 


Wing area, aq ft 


1, 626 


1,826 


Span, ft 




3/|0 


Mean ge one trio chord, ft 


.13 •Oh 


13. 0i| 


Natural wing frequency, ope 


2.50 


1.1*3 


Slope of lift curve, per radian 


1*.95 




Forward velocity, mph 


190 


160 






1*.97U.2 


v * 


60,882.6 


97,025.8 


Spring constant, lb per ft 


12.i405.75 


12.1405.75 


An 


2 


2 


Model B 








Condition I 


Weight, lb 




100,000 


Wing area, eq ft 




1.710 


Span, ft 




II4O 


Mb an geometric chord, ft 




12.21 


Natural wing frequency, ope 




2.1i5 


Slope of lift curve, per radian 




5.01; 


Forward velocity, mph 




260 


w e 




3.i»25.1* 



w fe , lb 96.57U.6 
Spring constant, lb per ft 85*233.1 
An 2 



TABLE II,- MOST SERIOUS VALUES OF OVERSTRESS FROM ADDING 
THE REACTIONS OF TWO GUSTS 



Hi 


H 3 


H 3 


Maximum 


chord 


chord 


chord 






lengths 


lengths 


SdAst 




Model A, condition I 


1+.26 


1+.26 


H+.50 


1 On 


1+.26 


10.25 


11.96 


l.OQ 

4- # w y 


1+.26 


19.53 


5.13 


1.09 


10.25 


I+.26 


35.01+ 


T 71 


*10.25 


10,25 


32.1+8 


1.21 


10.25 


19.53 


25.61+ 


1.20 


19.53 


1+.26 


69.22 


l.i+2 


19.53 


10.25 


66.66 


1.31+ 


19.53 


19.53 


59.82 


1.39 


Model A, condition II 




3.96 


3.96 


15-1+7 


1.1+7 


3.96 


10.08 


13.31 


1.35 


3.96 


20.15 


11.87 


1.35 


10.08 


3.96 


35.98 


1.1+7 


*10.08 


10.08 


33.83 


1.26 


10.08 


20.15 


32.39 


l.ll 


20.15 


3.96 


73.1+1 


1.56 


20.15 


10.08 


71.25 


1.36 


20.15 


20.15 


69.8I 


1.23 


Model B, condition I 


3.75 


3.75 


16.23 


1.1+0 


3.75 


9.99 


13.71+ 


1.26 


3.75 


19.98 


3.75 


1.26 


9.99 


3.75 


39.91+ 


1.35 


*9.99 


9.99 


57.1+6 


1.25 


9.99 


19.98 


27.1+7 


1.08 


19.98 


3.75 


72.1+3 


1.1+1+ 


19.98 


9.99 


69.93 


1.31+ 


19.98 


19.98 


59.91+ 


1.20 



*Signiflcaiice given in text. 




In tht calculations C *3 dis regarded <*nd } in 
light condition j^w^ff =-0 

Figure I.- Equivalent wtng aM spriflQ system* 
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Fig. 4 




2.0 2.4 2.6 3.2 3.6 
. . „ Time, sec 

(a) Hi = 3.96 chordsjb = 3.91. (b) Hx = 10.08 chords;b = 1.30. (c) Hi = 20.15 chordsjb = 0.494. 

Figure 4.- History of reactions in a single gust. Model A, condition II. 
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Figs. 6,7 




0 -4 .8 1.2 1.6 20 2.4 2.8 3.2 3.6 4.0 

Time, sec 

Fi ur^ •">.- Fistor-' of reactions in'a repeated ^u?t. ?.'odel A, condition I, 
V-l = H 2 = 10.25 chord?, V, = 32.4* chord s,U 2 = -t^. 
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0 4 8 12 16 20 

H t ,chord lengths 

Curve »,D J $ W /An r Curve B, 6 d /& st Curve Ct&WAn- 

Figure 9.- Variation of ratios of maximum dynamic to static 

reactions with gradient distance, Model A, condition I. 
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4 8 12 16 20 

H 1 .chord lengths 
Curve \,D 2 S W /An r Curve B, 6 d /6 8t Curve 0,D*<5f/An P 

Figure 10.- Variation of ratios of maximum dynamic to static ? 

reactions with gradient distance. Model A, condition II £ 




.60 



Computed points 
-Model A ^condition I 
- "A, « If 
« B t " / 




.65 JO 
Ratio ofU rpi to U mo 



Figure 13.- Variation of effective stress ratio with ratio of 
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Curve A,0 S w /An r Curve B, 6 d /6 8t Curve C,D 2 S f /An r 
Figure 11.- Variation of ratios of maximum dynamic to static 
- reactions with gradient distance. Model B, condition I. 



